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Earlier  predictions  of  toughness  by  Thomson  and  Weertman  are  reviewed. 

The  differing  predictions  of  the  two  authors  are  shown  to  be  due  to  different 
interpretations  of  the  functional  dependence  of  the  size  of  the  elastic  region 
on  intrinsic  surface  energy.  An  analysis  of  the  quasi  continuum  model  in  terms 
of  dislocation  rearrangements  Is  given,  and  the  nature  of  the  boundary  condi¬ 
tion  on  the  elastic  enclave  boundary  is  discussed. 
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I.  Introduction 


Thomson1  and  Weertman  have  recently  published  expressions  (1)  and 
(2),  respectively,  for  the  stress  intensity  factors  for  equilibrium  cracks 
based  on  an  elastic  enclave  model  of  fracture 
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(We  have  modified  Eqn.  (1)  from  that  given  in  Ref.  1  by  a  numerical  factor 
of  2^1+n^Zn  to  conform  to  the  analysis  carried  out  below.) 

K  is  the  stress  intensity  factor  as  measured  at  the  external  grips  for  a  crack 
in  small  scale  yielding  of  a  given  length,  o  the  yield  strength  of  the 
material,  p  the  shear  modulus,  r  the  intrinsic  surface  energy,  n  the  strain 
hardening  parameter,  R  the  radius  of  the  elastic  enclave,  v  Poisson's  ratio, 
and  c  a  derived  quantity,  c<<1.  They  differ  most  strikingly  in  that  Eqn.  (1) 
predicts  a  power  law  dependence  on  r  with  a  power  in  the  neighborhood  of  2-3, 
while  Eqn.  (2)  apparently  yields  the_linear  dependence  of  pure  brittle  fracture, 
except  for  the  multiplying  factor,  c~2.  We  will  show  that  one  may  find  a  range 
of  dependence  on  r  from  one  of  these  expressions  to  the  other,  depending  on  how 
one  interprets  the  variation  of  the  inner  radius  R  with  external  stress.  The 
two  expressions  were  derived  using  different  boundary  conditions  on  R  ,  however, 
and  a  second  aim  of  the  paper  will  be  to  explore  qualitatively  the  question  of 
what  the  boundary  conditions  should  be  in  terms  of  dislocation  models. 

We  will  work  exclusively  in  Mode  III  anti-plane  strain  for  the  reason  that 
the  mathematical  analysis  in  this  case  is  simple,  and  the  consequences  for  a 
specific  physical  model  can  be  clearly  worked  out.  We  assume,  of  course,  that 
predictions  for  Mode  III  should  have  qualitative  significance  for  the  physi¬ 
cally  important  Mode  I  case.  Also,  we  shall  work  with  a  time  independent, 
static,  deformation  law  for  the  reason  again  that  the  crucial  issue  is  to 
decide  how  to  model  the  crack  physically  with  an  elastic  enclave,  and  not  to 
make  detailed  predictions  regarding  crack  growth. 

When  we  think  of  embedding  a  sharp  crack  In  a  material  containing  external 
dislocation  sources  distributed  randomly  around  the  crack  tip,  the  tip  of  the 
crack,  on  probabilistic  grounds,  is  most  likely  not  to  experience  the  blunting 
catastrophe  which  would  occur  If  a  slip  plane  should  intersect  the  exact  position 
of  the  tip  over  a  significant  front  of  the  crack  and  destroy  its  sharp  character. 
Rather,  slip  taking  place  ahead  of  the  crack  will  pass  it  by  without  changing 
the  geometric  configuration  of  the  cohesive  zone,  and  slip  taking  place  behind 
it  will  cut  the  cleavage  surface  and  create  typical  slip  steps  on  the  cleavage 
surface.  The  crack  will  then  be  supposed  to  take  the  general  shape  depicted 
In  Figure  1.  It  is  our  object  to  work  out  the  predicted  toughness  of  a 
material  possessing  such  a  crack. 

II.  The  Elastic  Enclave  Model  in  Mode  III 


To  derive  any  quantitative  results,  it  Is  necessary  to  model  the  plasti¬ 
cally  deformed  region  surrounding  the  crack  by  a  continuum  theory.  Further, 
in  the  usual  fashion,  we  assume  small  scale  yielding,  and  assert  that  outside 
some  finite  boundary  surrounding  the  crack  the  material  Is  stressed  below  its 
yield  stress,  and  is  linear  elastic.  We  label  these  regions  In  Figure  2,  II 
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Figure  1  -  A  sharp  crack  which  has  been  blunted  by 
impingement  of  discrete  slip  lines  on  its  cleavage 


m  ELASTIC  CONTINUUM 

Figure  2  -  Schematic  of  a  crack  surrounded  by  elastic 
and  plastic  regions.  I,  elastic  enclave  surrounding  the 
crack  tip  of  radius  R  ;  II,  plastic  continuum  region  of 
radius  R  ;  III,  elastic  region.  The  crack  Is  centered 

In  R_,  but  R  Is  not  concentric, 
e  p 

and  III,  respectively.  Finally,  If  the  Immediate  local  vicinity  of  the  crack 
tip  possesses  no  dislocation  sources  and  Is  devoid  of  slip  planes,  then  In 
this  local  region,  termed  by  Weertman  the  elastic  enclave,  the  material  Is 
again  linearly  elastic  until  the  stresses  become  so  high  In  the  cohesive 
region  of  the  crack  that  discrete  nonlinear  effects  must  be  considered.  In 
the  following,  we  think  of  the  cohesive  region  as  characterized  by  a  critical 
local  K  in  the  Barenblatt  sense,  and  consider  the  crack  to  be  an  elastic  singu¬ 
larity  with  this  K.  In  Figure  2  we  label  the  elastic  enclave  I. 


51 


For  reasons  which  will  become  apparent  in  the  following,  we  have  shown 
all  boundaries  as  circles.  The  outer  boundary  between  II  and  III  is 
characterizec  as  the  locus  of  points  where  the  stress  equals  the  critical 
shear  stress,  a  .  The  inner  boundary,  however,  is  more  difficult  to  specify. 
In  reference  (l“  we  related  it  to  the  maximum  dislocation  density.  However, 
the  meaning  of  region  I  is  best  seen  by  considering  the  qualitative  form  of 
the  crack  opening  displacement  as  a  function  of  the  distance  from  the  crack 
tip.  In  Figure  3  we  show  in  a)  the  displacement  function  for  the  cleavage 
surface  for  a  fully  brittle  continuum  crack,  and  in  b)  that  for  a  continuum 
plastic  crack.  In  c)  we  show  what  a  crack  displacement  function  would  look 
like  if  the  deformation  is  heterogeneous.  From  these  figures,  we  should 
identify  R  with  the  distance  between  the  discrete  steps  representing  actual 
slip  planes  near  the  crack  tip.  Further,  as  a  plausible  procedure  for  solving 
the  stress  distribution  problem,  we  shall  interpret  R  as  a  lower  "cut-off" 
for  the  plastic  solution,  and  assume  that  the  solution  for  R  >  R  is  given  by 
the  continuum  plastic  solution.  Within  Rg,  only  the  elastic  portions  of  the 
strain  and  displacement  functions  persist!  and  we  match  elastic  values  of  the 
solution  on  the  boundary  R  .  Since  there  are  no  external  sources  of  stress 
on  R  ,  the  stress  must  also  be  continuous  there. 


Figure  3  -  a.  crack  opening  of  elastic  crack;  b.  crack 
opening  of  fully  plastic  crack;  c.  crack  opening  when 
the  plastic  strain  is  discretized  on  slip  planes. 

Recapitulating,  our  procedure  amounts  to  a  recognition  that  the  displace¬ 
ment  Is  composed  of  two  components,  elastic  and  plastic,  which  normally  must 
be  coupled  to  determine  the  overall  solution,  but  in  our  case  must  be  decoupled  • 
because  the  plastic  component  undergoes  a  cut-off  before  the  elastic  component 
does.  Within  the  confines  of  a  continuum  theory,  there  is  no  rigorous  way  to 
accomplish  this  decomposition,  but  clearly  an  approximate  solution  will  be 
achieved  If  we  adopt  the  full  continuum  solution  of  Figure  3b.,  separate  off 
the  elastic  portion,  and  match  an  elastic  crack  in  I  to  the  elastic  component 
of  II  at  Rfi. 

3 

Analytically,  we  follow  Rice's  solution  of  the  work  hardening  Mode  III 
problem.  We  assume  a  constitutive  work  hardening  law  of  the  form 

0  ‘  (3) 

To 
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where  a  and  y  are  principal  stress  and  strain,  respectively,  and  a  ,  and  y 
are  the  corresponding  critical  yield  values.  We  assume  the  crack  Ties  along 
the  negative  x-axIs  with  Its  tip  at  the  origin.  In  Rice's  problem  R  =  0.  He 
Introduces  a  potential  function  ,<j> ) .  where  $  is  the  angle  between  y  and  the 
y  axis,  with  the  defining  equations. 

X  =  -sin  4>  3V  -  COS  t  3V 

dy  y  3-t 

w  =  cos  *  3f  -  sin  *  3¥ 

¥7  y  H  (4> 

which  satisfies  the  differential  equation 


Y2  d2*  +  ny  11  +  n32V  _ 
3y*  3y  3$  2 

The  boundary  conditions  are  that 


III 


ip  =  ±n/2 


=  0 


(5) 


(6) 


If  we  solve  Eqn.  (5)  by  separation  of  variables  then  v  must  have  the  form 

(7) 


v  ■  v(y)  sin  m  <f> 
m  =  ±1 ,  ±2— 


Then  Eqn.  (5)  is  solved  by  setting 

V  =  y*5 
m  ' 

where  p  satisfies  the  quadratic  Indiclal  equation 
p(p-l)  -  np  -  nm2  =  0 

For  each  value  of  m,  p  has  two  solutions,  p+  and  p_. 
The  general  solution  for  Eqn.  (4)  Is  the  series 


(8) 

(9) 


»  P+(m)  p  (m) 

T  =  t  (amY  +  bmr  )  sin  m  * 
m»l  ' ' 


(10) 


Rice  was  able  to  satisfy  the  boundary  conditions  for  iu  >  0  at  the  crack 
tip  and  the  boundary  conditions  on  the  elesi.1c-plast1c  boundary  by  retaining 
only  one  term  of  the  series  in  Eqn.  (10), 

»  *  a,  sin  ♦  (11) 
yn 


From  this,  he  shows 


x  ■  X(y)  +  R(y)  cos  2  * 
y  *  R(y)  sin  2  p 
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02) 


X(y)  =  ^  R(y) 


2$  =  e 


This  solution  is  represented  by  a  set  of  nested  circles.  On  the  circumference 
of  each  circle,  y  is  a  vector  of  constant  magnitude,  whose  polar  components  are 
given  by  the  equations 


yR  =  Y  sin  e/2  oR  =  a  sin  e/2 
Y„  =  Y  cos  9/2  a.  =  a  cos  9/2 
9=2$ 


03) 


where  R  and  e  are  measured  in  a  system  of  polar  coordinates  with  the  origin 
at  the  center  of  each  circle.  For  different  values  of  y,  the  centers  of  the 
circles  shift,  being  situated  at  x  =  X(y),  relative  to  the  crack  tip.  See 
Figure  4.  The  stress  intensity  factor,  K,  defines  the  radius  of  the  elastic 
plastic  boundary,  Rp 

K2  =  2"a02Rp  (14) 


Figure  4  -  Schematic  drawing  of  crack  embedded  in  three 
regions  as  in  Figure  2.  The  accumulation  point  for  the 
plastic  circles,  R(y),  Is  to  the  left  of  the  center  of 
the  elastic  enclave.  X (y )  is  the  distance  from  the 
accumulation  point  to  the  center  of  the  circle,  R(y). 

An  arbitrary  R  (y)1s  shown  with  the  principal  strain 
y  on  the  circumference,  and  coordinates,  R,  9. 
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Figure  5  -  The  continuum  plastic  solution  in  our  model 
is  replaced  by  an  elastic  crack  opening  at  the  tip  with 
a  discontinuous  displacement  ar  R  corresponding  to 
concentrating  all  the  continuum  dislocations  within  R 
onto  its  circumference.  The  slope  of  the  displacement 
at  R  is  continuous. 


Equations  (4)  through  (14)  describe  a  fully  continuum  plastic  crack.  In 
our  problem,  we  assume  at  a  value  of  R  given  by  R  ,  the  plastic  solution  is 
cut  off,  and  the  solution  is  continued  within  R  By  a  purely  elastic  solution. 
(See  Figure  5.)  We  now  write  the  solution  for  a  purely  elastic  crack  situated 
at  the  geometric  center  of  R  .  K’th  the  usual  boundary  conditions  au  =  0  at 
the  tip,  and  zero  stress  on  the  crack  surfaces,  the  solution  within  R  is 


l 

c=2 


n/2 


z  *  x-iy 

o(z)  =  o32(x,y)  +  i  o31(x,y)  (15) 

where  Cn  is  real,  and  K  is  the  local  K  for  the  crack  tip  enclosed  in  the 
elasticnconclave.  In  polar  coordinates  centered  at  the  center  of  R  ,  the 
first  term  has  the  form  e 


-  -  SIM  -  UTn 

K  K 

^  R  <  R_ 

0Q  =  —  cos  e/2  =  UT. 

9  J7JH  8 

u  =  ^  sin  e/2  (16) 

where  u  is  the  component  of  the  displacement  normal  to  the  x-y  plane. 
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In  the  foregoing,  we  have  used  a  variety  of  coordinate  origins,  and  crack 
tip  positions,  which  we  wish  to  clarify  here  to  avoid  possible  confusion  for 
the  reader.  In  writing  the  equations  (4)  through  (14)  describing  the  plastic 
solution,  we  had  reference  to  an  effective  crack  tip  position  shown  in  Figure  4 
to  which  the  continuum  solutions  were  referred.  This  effective  tip  position 
was  the  accumulation  point  about  which  the  nested  circles  converged  as  R  ■*  0. 
Also,  we  described  a  coordinate  system  for  e  in  Eqn.  (13)  which  specified  the 
orientation  of  the  principal  strain  (and  stress)  for  a  particular  magnitude 
of  y  and  for  a  given  R,  where  R  >  R  >  R  .  This  coordinate  system  was  centered 
at  the  center  of  the  given  circle  R.  This  center  is  at  a  distance  X(y)  from 
the  effective  position  of  the  continuum  plastic  crack.  Finally,  we  assume 
that  a  crack  tip,  elastic  this  time,  is  placed  at  the  center  of  the  circle, 

R  ,  whose  circumference  defines  the  boundary  of  the  elastic  enclave.  We  note 
that  by  construction,  the  accumulation  point  of  the  nested  circles  of  the 
plastic  solution  at  the  effective  plastic  crack  position  is  not  at  the  same 
position  at  the  elastic  crack  at  the  Center  of  R  .  Nevertheless,  the  actual 
physical  crack  is  assumed  to  be  positioned  at  the  center  of  R  ,  and  coincides 
with  the  elastic  crack  of  equations  (15)  and  (16).  Its  solution  within  R  is 
given  by  equations  (15)  and  (16),  and  outside  R  is  given  by  equations  (4?-(14) 
with  boundary  conditions  given  on  Rfi  which  we  now  specify. 

On  the  circle,  R  ,  we  set  the  stress  of  the  elastic  solution,  Eqn.  (16), 
equal  to  the  stress  o“  the  plastic  solution,  Eqn.  (13).  We  note  that  since 
these  equations  display  the  same  angular  dependence,  with  respect  to  e  on  R  , 
this  is  easily  done.  With  the  use  of  the  third  equation  of  Eqn.  (12)  and  tfie 
defining  relation  for  K,  we  have 


(17) 


If  we  assume  that  the  elastic  crack  is  in  equilibrium,  that  is,  K  obeys 
the  Griffith  relation,  e  J 


Ke2  =  4  ru 


we  have 


R  (l-n)/(1+n) 


K2  =  4  ru  (A 

Ke 

*  2*  (^)(1  ",/(Zn)  (Reo02)  (1'n)/(2n> 


(18) 

(19a) 

(19b) 


This  Is  the  result  quoted  In  Eqn.  (1). 

*5**  :esu1t  "eertman's  results,  Eqn.  (2),  which  we 
take  from  Eqn.  (24)  of  reference  (2).  If  we  rewrite  this  equation  In  our 
notation,  we  have 


a  (l-n)/(2n) 

C  "  n  (5TO) 


(20) 


Weertman’s  parameter  p  can  be  replaced  by  1/(n+l),  and  n  is  a  constant,  n<l. 
Since  in  Node  III,  (o0MRe))  -  (Rp/Re)n/n+1,  we  have 
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(21) 


K  (l-n)/2(l+n) 
c  =  n  (A 
RP 

Except  for  the  constant,  n,  this  is  precisely  the  result  (19a)  of  this  paper. 

Thus,  a  direct  comparison  of  Eqn.  (1)  with  Eqn.  (2)  shows  they  are  based 
on  precisely  the  same  results.  They  differ  in  that  in  Eqn.  (1),  we  have  not 
taken  the  ratio,  R  /R  ,  to  he  a  constant.  There,  R  is  written  in  terms  of 
its  dependence  on  r,  ind  R  is  left  as  an  undetermined  parameter.  Thus  the 
discrepancy  in  the  published  results  apparently  does  not  revolve  around  a 
different  use  of  the  boundary  conditions,  but  rather  depends  upon  how  the 
ratio  Rp/Re  is  interpreted. 

In  our  opinion,  there  is  no  physical  reason  why  R  and  R  will  depend  on 
r  in  precisely  the  same  manner.  Hence,  the  dependencepof  K  ofi  r  should  be 
more  complex  than  a  simple  linear  relation.  From  the  discussion  here,  however, 
it  is  clear  that  no  theory  yet  exists  from  which  R  can  be  related  to  r,  and  in 
fact,  from  Eqn.  (19a),  experimental  measurement  ofeK(r)  would  provide  evidence 
on  the  point. 


III.  Alternative  Continuum  Solutions 


In  this  section,  we  explore  Weertman's  suggestion  that  the  total  dis¬ 
placement  at  R  should  be  continuous,  as  well  as  o.  However,  rather  than 
follow  Weertman's  Mode  I  development,  we  believe  it  is  preferable  to  remain 
in  Mode  III,  where  the  analysis  is  more  accessible,  and  compare  Weertman's 
prescription  with  the  results  of  Section  II  of  this  paper. 

Following  Weertman,  we  write  a  solution  in  region  II  of  the  form 

'rII“R~P,  and  then  match  appropriate  power  series  solutions  in  regions  I  and 
III  to  this  region  II  plastic  solution.  We  believe  that  a  simplification  is 
possible  which  is  suggested  by  the  fact  that  since  RD>>Re>  the  solutions  at 
R  and  the  matching  process  there  should  be  effectively  independent  of  the 
details  of  the  matching  on  R  .  That  is,  microstructural  effects  are  fully 
damped  out  when  R  ^  R  .  Thus,  in  Region  III,  we  follow  Rice  ,  and  take  the 
elastic  solution  in  III  to  be 
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III 
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—  sin  0/2 

SEX 
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II! 
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—  cos  e/2 


R  >  Rp  (22) 

R  and  a  are  the  polar  coordinates  relative  to  an  origin  fixed  at  the  center 
of  R  .  That  is.  In  region  III,  the  crack  tip  has  an  effective  position  at 
the  center  of  the  plastic  circle,  R  .  Also,  the  plastic  solution  at  R  should 
consist  only  of  the  expressions  fort  used  by  Rice  ,  which  was  the  first  term 
m=l,  of  the  expansion,  Eqns.  (7)  and  (10). 


However,  if  we  are  to  match  the  total  displacement  and  stress  at  R  ,  we 
shall  have  to  consider  additional  terms  In  either  the  expansion,  Eqn.  (75) , 
for  the  elastic  region  I  or  for  the  plastic  region  II,  or  both.  Further,  If 
we  Include  higher  order  terms  in  Eqn.  (10),  they  must  disappear  with  R  faster 
than  the  dominant  term,  „  R-l/(l+n).  An  inspection  of  the  two  expansions, 

Eqns.  (15)  and  (10),  however,  shows  that  in  order  to  match  solutions  on  some 
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boundary  corresponding  to  R  ,  a11_  the  terms  in  either  or  both  of  these  expan¬ 
sions  will  have  to  be  used,  because  their  angular  dependencies  are  inconmen- 
surate.  In  addition,  the  matching  surface  may  have  a  complex  shape  as  well. 
This  was  a  point  noted  by  Weertman,  who  in  the  face  of  this  difficulty,  simply 
ignored  the  angular  dependence,  and  matched  stress  and  displacement  with  the 
smallest  possible  number  of  terms  in  the  R  parts  only.  Following  this  path, 
we  can  again,  as  Weertman  did,  keep  only  the  dominant  solution  in  region  II, 
and  add  one  term  to  the  elastic  expansion  in  region  I.  at  R  =  R  we  have 


II  _  A,  "/<1+n>  a1  = 
0  -°o^ 


— —  +  A  /IT" 

vEHL  e 


II  Rn  ]/(1+n) 

u11  -  2y0  (A  Rfi 

e 


ui  =  A  +|  a 

U  *  71  J  U 


A  R 


3/2 


(23) 


These  expressions  are  obtained  from  Eqns.  (12),  (13),  C20),  and  (16),  but  we 
have  added  the  term  in  A  to  the  dominant  Kg  terms  in  <j  and  u  .  Writing 


&1-" 


+  A 


<£>'**  « 


/ +  i  a 

V  Tr  3  y 


„3/2 


(24) 


for  the  matching  conditions  on  Rfi,  we  finally  obtain 


1-n 

1+n 


1 
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Re 


1-n 

1+n 


(25) 


Since  R  »R  ,  Ecln-  (25)  shows  that  the  elastic  crack  is  not  shielded  by 
the  deformation  field,  and  K  >K,  which  Is  a  non-physical  result.  This  anomaly 
is  caused  by  the  large  value  of  the  elastic  K-field  required  to  match  the 
large  continuum  displacement  at  R  .  Also,  A  has  a  large  negative  value, 
which  requires  a  correspondingly  Targe  counterbalancing  K-field. 


IV.  Conclusion 


The  published  results  of  Thomson  and  Weertman  are  shown  to  reflect 
equivalent  physical  pictures  at  the  crack  tip.  Their  difference  is  caused  by 
different  implications  concerning  the  material  parameters  which  control  the 
size  of  the  elastic  enclave.  However,  as  there  has  yet  been  no  adequate 
treatment  of  the  size  of  the  elastic  enclave,  we  believe  It  would  be  fruitful 
at  this  stage  to  highlight  it  as  a  fruitful  field  for  experimental  study. 


In  effect,  what  our  treatment  does  is  to  modify  the  dislocation  distribution 
and  their  sources  In  the  region  of  the  crack  tip  from  that  of  the  plastic 
continuum  solution.  The  dislocation  distribution  Implied  by  the  plastic 
solution  Is  often  termed  the  "geometric"  distribution,  and  can  be  calculated 
from  the  curl  of  the  local  stress  function,  o.  The  geometric  distribution 
calculated  In  this  way  does  not  have  a  strong  singularity  at  the  crack  tip, 
and  It  Is  therefore  possible  to  modify  It  by  concentrating  the  distribution 
totally  on  the  circumference  of  R  without  markedly  changing  the  overall  K. 
However,  modifying  the  manner  in  Which  these  dislocations  are  Introduced  Into 
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the  medium  from  the  cleavage  surface,  for  example,  by  assuming  that  real 
dislocations  form  discrete  slip  lines  on  crystallographic  planes,  will  modify 
the  detailed  shape  of  the  crack  opening  near  its  tip.  This  is  precisely  the 
point  of  our  discussion  as  shown  in  Figures  4  and  5,  of  course.  On  the  other 
hand,  these  modifications  in  crack  opening  near  the  tip  do  not  change  the 
total  crack  opening  displacement  (COD),  provided  the  total  integrated 
burgers  vector  of  the  screening  dislocation  distribution  is  unchanged,  as 
shown  again  in  Figures  5  and  6.  Finally,  such  changes  within  R  amount  to 
a  change  in  the  local  stress-strain  law,  which  will  modify  the  stress  and 
strain  solutions  within  R  markedly.  Indeed,  in  our  solution,  we  have 
recognized  this  importantchange  by  specifying  an  elastic  stress-strain 
law  within  R  ,  with  the  results  derived  in  section  II.  A  more  adequate 
theory,  of  course,  would  deal  in  a  unified  way  with  the  elastic  and  plastic 
regions  by  means  of  actual  discrete  distributions  of  dislocations.  Thus, 
although  we  believe  we  have  some  of  the  physics  right  in  Eqns.  (19a  and  19b), 
we  are  less  sure  of  the  details,  and  we  have  lumped  some  large  questions  into 
our  parameter,  Rg. 

In  summary,  we  do  not  believe  we  yet  have  an  adequate  theory  of  the 
toughness  exhibited  by  sharp  cracks  in  a  medium  of  modest  ductility.  The 
problem  remaining  relates  to  how  one  makes  a  rigorous  transition  from  the 
plastic  continuum  to  the  discrete  description  of  dislocation  distributions 
at  the  crack  tip.  The  cut-off  procedure  of  this  paper  is  at  best  ad  hoc,  and 
since  the  cutoff  radius,  R  ,  appears  in  an  essential  way  in  the  results  (i.e., 
leaving  open  the  dependence  of  R  on  r,  etc.),  we  believe  a  proper  theoretical 
treatment  still  lies  in  the  future. 

We  have  profited  from  discussions  of  this  problem  with  E.  Hart,  I.  Lin, 
and  J.  Weertman. 
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